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Algorithms that reject transient signals due to proton effects on charge coupled device (CCD) sensors have 
been implemented in the HDOS ASTRA-1 Star Trackers to be flown on the TOPEX mission scheduled for 
launch in July 1992. A unique technique for simulating a proton-rich environment to test trackers is described, 
as well as the test results obtained. 
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Solar flares or an oibit that passes through the South Atlantic Anomaly can subject the vehicle to very high 
proton flux levels. There are three ways in which spurious proton generated signals can impact tracker per- 
formance: the many false signals can prevent or extend the time to acquire a star; a proton-generated signal 
can compromise the accuracy of the star’s reported magnitude and position; and the tracked star can be lost, 
requiring reacquisition. Tests simulating a proton-rich environment were performed on two ASTRA- 1 Star 
Trackers utilizing these new algorithms. There were no false acquisitions, no lost stars, and a significant re- 
duction in reported position errors due to these improvements. 


STAR TRACKER OPERATION IN A HIGH DENSITY PROTON FIELD 


Star sensors have been utilized in high-accuracy attitude determination systems since the early 1960’s, pro- 
viding precise measurement of the position and magnitude of stars in the sensor’s field of view (FOV). As 
applications continue to demand higher performance, the effects of a natural or enhanced radiation environment 
need to be accommodated by the star sensor’s design. HDOS recently delivered two star trackers for the 
TOPEX/Poseidon mission which is scheduled for launch in July 1992. These enhanced trackers contain an effi- 
cient mix of hardware and firmware that permits effective acquisition and tracking throughout their orbit, which 
includes extensive exposure to the transient rich environment of the South Atlantic Anomaly. This paper dis- 
cusses the algorithms employed, the environment simulated, and the results of the tests performed, demon- 
strating successful operation in a proton-rich environment. 


BACKGROUND 

State-of-the-art star sensors (see Figure 1) have recently benefited from two major technological devel- 
opments, the CCD detector and the microprocessor. The heart of the current HDOS ASTRA star sensors is 
the RCA 504 CCD a 256 x 403 pixel array which operates in the frame transfer mode. The thinned, backside 
illuminated device provides high quantum efficiency in the visible range. A thermoelectric cooler is used in the 
ASTRA-1 sensors to provide low noise operation in harsh environments. Fitted with a wide FOV (7x9 de- 
cree) color corrected lens, the ASTRA star sensors can provide position accuracy to 10 arc-seconds or better, 
and sensitivity down to a visual magnitude of 6. The HDOS ASTRA star sensors also utilize a versatile 16- 
bit microprocessor. Acquisition and tracking, centroid determination and correction, debris and transient event 
discrimination, and self-test functions are performed autonomously by the microprocessor. This design pro- 
vides a flexible interface and reduces the computation burden on the host computer. 

Star sensors utilizing mosaic CCD arrays can be separated into two classes: star trackers and star mappers. 
Star trackers have two distinct functional modes: acquisition and track. During acquisition, the sensor field of 
view is searched for valid targets. High-pass spatial filtering and pixel amplitude thresholding are often used 
to limit the amount of data saved each frame. Data from multiple frames can be used during acquisition to 
discriminate valid stars from debris and transient events. Once a star has been acquired, the sensor enters 
the track mode. During the track mode, data from previous frames is used to estimate the current position of 
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Figure 1. Functional Diagram of the Advanced Star Tracker (ASTRA). 

the star. A reduced field surrounding the star is then read, processed, and the magnitude and position of the 
star are output to the host. The HDOS ASTRA- 1 built for the TOPEX/Poseidon mission is an example of a 
star tracker. The ASTRA Star Tracker currently being built for Space Station Freedom will operate as four 
virtual star trackers allowing simultaneous acquisition and tracking of up to four stars. 

Star mappers acquire stars throughout the total sensor field of view and report their magnitude and position 
each frame. High pass filtering and pixel amplitude thresholding are used to limit the amount of data that must 
be processed each frame. Normally, information is not saved from frame to frame. The HDOS ASTRA-2 built 
for the JH/APL START experiment 1 is an example of a star mapper. It outputs the position and magnitude of 
up to five stars at a 10 Hz update rate. 


Trackers have a number of functional advantages over star mappers when operating in a transient event-rich 
environment. Once a star tracker has acquired a star, only a small fixed number of pixels must be read, stored, 
and processed each frame. A star mapper must process the entire CCD array each frame and the number of 
transient events will impact the amount of data that must be processed. By using data from multiple frames, 
star trackers can discriminate between stars (with predictable magnitudes and positions) and transient 
events during acquisition. Using magnitude and position data from previous frames during track mode also al- 
lows a star tracker to determine if a transient event has corrupted the reported magnitude and position data. 
Star mappers do not store information from previous frames and cannot discriminate between transient events 
and valid targets; this task must be performed by the host. Star mappers can also report a different set of 
stars each frame which requires identification of stars each frame. In a star tracker system, once a valid star 
has been identified, there is no need for the host to re-identify a star so long as it remains in track. 

THE ENVIRONMENT 

For a star sensor to operate in a natural or enhanced radiation environment, the effects of radiation damage 
and radiation induced noise events on the CCD must be addressed. Space radiation that will interact with the 
CCD can be divided into two major groups, trapped particles and solar protons. Trapped particles are pro- 
tons and electrons trapped in the magnetic fields surrounding the earth. The energy and fluence profiles will 
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vary with the altitude and inclination of the orbit. In low earth orbits, anomalies at the poles and off the coast 
of Brazil (the South Atlantic Anomaly) can result in high fluxes of energetic protons. Solar protons are emitted 
during solar flares which occur randomly. Intense solar flares occur about every 11 years. At high altitudes, 
i.e. geostationary orbit, solar protons can dominate, with fluxes sometimes orders of magnitude greater than 
the trapped particle radiation. HDOS has developed the capability to evaluate the radiation 10 

space and its impact on sensor systems for given orbits. These models are based on the RADBELT AP-8, 
AE-8 radiation environment data supplied by NASA-GSFC. 

The effects of radiation damage in CCDs have been the focus of much recent work. 4 HDOS has developed the 
tools required to model the end-of-life performance of CCD detectors, and has also performed experiments to 
measure these effects on irradiated devices. The primary concern of this paper is the effect of transient noise 
events that occur when high energy particles impact the CCD. Charged particles interact directly with the 
CCD pixels, causing a random series of ionization events. These events can be localized to within a few pix- 
els or can result in streaks, depending on the angle of incidence of the particle, the geometry of the device, and 
the energy of the particle. 5 6 A large number of events in the sensor field of view can interfere with the opera- 
tion of the sensor. If the system must operate in a high density proton environment, a method to reduce the 
impact of these events on the sensor operation is required. The most obvious solution is to increase the 
shielding around the CCD to reduce the number of events. However, this will increase system size and 
weight. A more elegant solution is to apply real-time processing to reject these transient events. 

Transient events can degrade a star tracker’s performance in a number of ways. During acquisition, transient 
events can be acquired falsely or can impede acquisition of a valid target. During track, transient events can 
corrupt position and magnitude data or can result in the sensor dropping a valid star track. For the TOPE 
mission, P Fairchild Space Co. 7 determined that the star tracker must operate with up to 150 transient events in 
the sensor FOV per frame, 100 milliseconds. All of these events were assumed to be indistinguishable from 
real stars to the CCD. Since these transient events were independent of one another their position and 
magnitude were random and uniformly distributed over the CCD array. 


The system had to meet the following requirements: 

• Acquire and track stars with up to 150 transient events per frame 

• Acquire and track stars moving up to 0.3 degree/second 

• Probability of acquiring a valid star within 22 seconds is 95% 

• Alert host if data has been corrupted 

• Maintain track of valid stars during proton event disturbances. 


ALGORITHM IMPLEMENTATION 

Any solution is dependent on determining what information is available to discriminate between transient 
events and valid stars. Star image size is a function of the point spread of the optics and the motion and jitter 
of the spacecraft. Therefore, bounds can be placed on the size of valid targets and single pixel upsets and 
events that result in long streaks can be rejected. Proton events only last for one frame. Since they are in- 
dependent of one another, the position and magnitudes of the proton events randomly change from frame to 
frame Positions of the star images on the array are determined by the dynamics of the spacecraft and there- 
fore change systematically from frame to frame. By applying spatial and magnitude filtering, we can reject 
transient events. The probability of transient events passing this filter will be a function of the number o 
transient events per frame and the size of the discrimination windows. The line-of-sight (LOS) motion re- 
quirement (up to 0.3 degree/second) thus determines the minimum spatial window to be used. The size ot the 
magnitude window is determined by the tracker’s predicted error in determining star magnitude. Once a star 
is acquired, magnitude can be checked to determine if a transient event is contiguous with the star image, 
causing an error in the centroid of the image and in its reported magnitude. 
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Our selected approach was to limit the data processing load during each tracker frame when in the acquisition 
mode to be consistent with existing hardware capabilities. The tracker FOV was divided into 19 acquisition 
bands (19 pixel rows by 403 pixel columns). The bands were stacked in the row direction with a small overlap 
and the search for a star confined to one band at a time. If the acquisition sequence fails to find a star in the 
current band, the tracker sequences to the next band and searches for a valid star. This sequence continues 
until the tracker acquires a star in one of the bands and transitions to track mode. If the entire FOV is 
searched unsuccessfully, the sequence is repeated. 

The acquisition sequence consists of two acquisition states (A1 and A2) and five validation states (HI 
through H5). Each state takes one tracker frame to implement. 

The acquisition sequence searches an acquisition band to determine if any group of pixels with signal exceed- 
ing a threshold exists that could quality as a star. If no such group exists the tracker transfers to the next 
acquisition band and the sequence is repeated. If candidate stars are found in the A1 state, the tracker enters 
the A2 state. If any of the candidate stars identified in A1 is also found in A2, the tracker goes into the vali- 
dation phase of the acquisition sequence. During validation states, the tracker continues to evaluate the can- 
didate star’s characteristics for temporal and spatial consistency. Upon transition to the track state, sufficient 
history has been established so that it is highly probable that the group represents a valid star and is not 
caused by transient events. If the star is not confirmed in any of the validation states the tracker increments 
to the next acquisition band and the acquisition sequence is continued. 

Upon transition to track, a small track window is defined which is centered about the last validated star po- 
sition. The window position is updated each frame to track the updated star position. The window is made 
small enough to limit the data processing load and exclude the majority of transient proton signals; but it is 
large enough to account for vehicle LOS rates. 

The track window is scanned in a raster fashion each frame, and the magnitude and position characteristics of 
the candidate star are evaluated. If the star is evaluated as a “valid” star, the centroid data and star magni- 
tude data are sent to the host and a bit is set in the data word that indicates the data is valid. If no valid star 
is found in the track window after multiple attempts to recover from proton hits, the star is considered lost and 
the tracker reverts to the acquisition mode. 

TEST METHOD 

The CCD’s inability to distinquish between proton or photon generated-signals was utilized to test the 
tracker with the proton flux improvements. A Scene Simulator (SS) concept was implemented. Scenes con- 
sisting of point sources of light were generated on a computer monitor. The point sources were collimated and 
then imaged by the tracker on its CCD detector. The signals generated in this fashion at the CCD could be 
considered as having been generated by either protons or stars. This technique achieved complete control of 
the interactions between a simulated star and simulated proton events. 

The SS equipment shown in Figure 2 consisted of a VGA monochrome monitor and a personal computer (PC) 
with a 386 processor running at 25 MHz. The monitor was positioned at the focal plane of an achromatic colli- 
mating lens with a focal length (FL) of 1 185 mm. The scene was generated by commanding 150 monitor pixels 
on, to serve as “proton” sources, and one pixel to serve as a “star”. The intensity of each monitor pixel was 
controlled by the SS software. The collimated light from each monitor pixel was imaged by the tracker optics 
on its CCD detector. The tracker optics had an FL of 41 mm which imaged the monitor pixels at the CCD de- 
tector, resulting in a demagnification factor of 29. All images at the tracker, simulated protons as well as 
simulated stars, were equivalent in size to predictions for a valid star in orbit, i.e. between two and 25 CCD 
pixels. The angular FOV that the monitor accommodated exceeded the 9x7 degree FOV of the tracker. 

To ensure a worst-case test scenario, all proton events were sized (in pixels) within the range considered by 
the tracker to be acceptable stars. Figure 3 (frames 299 and 300) shows the excellent image size achieved 
with the Scene Simulator. 
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Figure 2. Configuration for Simulation of Stars and Proton Events in a 
Laboratory Environment. 
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Figure 3. Excellent Image Properties Were Achieved with the Scene Simulator (Mi~3.7). 

To limit simulated proton signals to one tracker frame, the generated scene display was s ^roni“d with die 
10 Hz frame rate of the tracker. This synchronization was mandatory because a star will be imaged at a spa 
ially cSLnt location during successive tracker frames, whereas a proton signal will be transitory both 

spatially and temporally. 


TESTS RESULTS 

The tests described were designed to verify three basic aspects of the tracker software performance in a pro- 
ton-rich field, specifically: 

• A star being tracked will not be lost as a result of the proton events 

• Acquisition of a valid star will occur within 22 seconds of the start of the search and star position 

• Magnitude data corrupted by proton events will be identified for the host. 
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The scene simulator technique was especially useful since the intensities and positions of the simulated pro- 
tons could be varied by the software to have any desired relationship to the simulated star. In all, 14 unique 
scenarios were used to test the response of the tracker to various combinations of “proton” influences on 
tracker performance. Generally the scenes simulated conditions of near or direct proton hits on the star posi- 
tion. In ^addition a no stars scene was used to verify that the tracker software was not fooled by the 
“proton” signals and did not erroneously report star acquisitions. Both moving and stationary stars were 
simulated. 

Each test scene scenario always consisted of 150 simulated protons and one simulated star. The simulated 
proton positions were generated in a random fashion in the tracker’s FOV. Consistent with Fairchiid’s radi- 
ation analysis/specifications, 35 of the simulated proton signals were at an equivalent star Instrument 
Magnitude (Mi) of 3.2; 115 ranged between the equivalent Mi of 3.7 and 5. The simulated star was set at 
Mi - 3.7. Generally the “star” in the scenes was not subjected to random “proton” influences. Each scene 
was carefully designed to introduce specific and periodic occurrences of “proton/star” interactions so that 
tracker software processes could be evaluated. Completely random occurrences of “proton/star” interactions 
were used to gather statistical information regarding frequency of “star” disturbances, to ensure that the 
three primary performance requirements were met and that interacting events did not cause unanticipated 
results. p 


The results of a test designed to demonstrate ASTRA’s capability to reject “protons” during the star acqui- 
sition sequence are shown in Figure 4. The particular scenario contains 150 “protons” but no valid star in the 
FOV. As required, at no time did the tracker indicate a star acquisition. If no signals were found during the 
search of the tracker’s FOV, the highest acquisition state achieved by the tracker is the A1 state and the 
minimum time to search the entire FOV is 3.8 seconds. The figure shows that the “protons” caused the 
search of each band to extend routinely to the A2 state, periodically to the HI state and occasionally to the 

H2 state. Each state beyond the A1 state, caused by proton signals, typically adds 100 msec to the overall 
search time. ' 


SceneOI O-Tst D1 220097.dat SceneOI O-Tst D1 220097.dat 



Figure 4. The ASTRA Design Never Reports a Proton as a Star. 

Figure 5 shows the results of a test in which a “star” moves horizontally across the FOV at an apparent ve- 
locity of 0.35 degree/second and periodically experiences a direct and “near proton” hit. The figure is similar 
to many of the following figures. Row, Column, and Intensity are displayed as a function of tracker frame 
number. Row and Column are plotted in CCD pixel space, Intensity is in output counts. Any time a star signal 
is invalid or lost, the values of all three parameters drop to zero. A momentary dropout is characterized by a 
zero signal for no more than 10 frames. If the signals go to zero for more than 10 frames, the tracker has 
reentered the acquisition mode, and the time that the signal remains at zero is indicative of the time required 
to reacquire the star. 
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Figure 5. ASTRA Does Not Report Star Data Corrupted by Direct or Very Near Proton Hits. 

Every 100 tracker frames, the SS causes a "proton” hit directly at the star position and 15 frames later a 
“Droton” hit occurs 0.035 degree from the star. It can be seen from the data in Figure 5 that each event 
causes a dropout (i.e., the star position data is not used by the host). The increased intensity of the direct 
“proton” hit causes the “star” to be evaluated as invalid since its increased intensity is outside the average 
intensity limits being continuously determined and updated by the tracker for the valid s ] , ar -,7^ 
frames later, drops out for the same reason since the near ‘ proton signal merges with the star si g 
increases the signal intensity. If the position determination had been reported for the near hit condition, the 
position would have been in error by the bias that the “proton” signal would have mtroduced. For the direct 
hit” condition an incorrect star magnitude would have been reported potentially impeding star identification by 

the host. 

It was also of interest to determine how close a “proton” could come to a “star” without perturbafing the 
star (see Figure 6). This scene has a diagonally moving star which is approached within 0.16 degree by a 
“proton” every 100 frames. As a benchmark, a simulated direct hit on the star by a proton was introduced 
causing a dropout for one frame, 15 frames prior to the near “proton” event. The data 15 frames after the 
dropout was reviewed and no effect of the close approach of the “proton” to the “star was evident. View (a) 
of Figure 6 shows the entire test results in which the regular “star” dropouts due to the direct P™ t0 
are evident. View (b) of the figure shows a single transit of the star over the tracker s FOV, view (c) shows 
a greatly magnified view of an area of interest. This plot is typical of all of the proton event occurrences dunng 

more than two hours of testing. 

Figure 7 demonstrates the ability of ASTRA to maintain track even when a significant number of sequential 
frames were impacted by interfering “protons”. The SS scenano was a stationary star interrup e y 
blank scenes every 100 frames. Despite the long periodic interruptions, the tracker did not revert to the ac- 

quisition mode. 

Table 1 provides statistical data on acquisition of a star in a “proton”-rich field. All the tests had either mov- 
ing stars that left the FOV at one edge and then were re-introduced at the other edge or stationary stars that 
were deliberately caused to break-track periodically. In all but the diagonal scan cases the SS scenano was 
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TRACKER PERFORMANCE COMPARISONS IN A PROTON FIELD 

ResuUs of tests in which the SS generates a stationary “star” positioned approximately in the center of the 
POV and 150 protons are generated in a completely random fashion in the FOV without regard to “star” 
position are shown in Figures 8 through 10. Random interferences with the “star” do occur. A comparison of 
the data from the three tests quantitatively demonstrates the increased error in the reported “star” position 
due to proton interference with the intensity discrimination disabled (Figure 9) and with both the intensity 
and position discrimination disabled (Figure 10). The latter case is indicative of the performance expected 
from the tracker if no consideration were made in the design to accommodate proton events. 
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Figure 7. Tracking of a Valid Star is Maintained for Multiple Sequential Proton Events. 

View (a) of Figure 8 illustrates the “star” position variation over an 11.6 minute test as reported by the 
tracker. An evaluation of the data shows that the uncorrupted data sent to the host represented 90 percent of 
the total number of frames. This statistic was valid for two separate trackers within 0.5 percent. This particu- 
lar test displayed a standard deviation of the reported row and column positions of 5.1 arc-seconds and 3.2 
arc-seconds, respectively. Figure 7(b) shows the row, column and intensity of the star plotted against time. 
Although there are numerous dropouts of data due to “proton” induced variations in star position or magni- 
tude beyond preset limits, the tracker never lost the “star” long enough for the tracker to re-enter the ac- 
quisition mode. 

Figure 9 presents the results from the same SS test scene but with the intensity compare discriminator dis- 
abled during the track mode. The standard deviation of the position data reported to the Host increased to 
14.1 and 12.5 arc-seconds in row and column positions, respectively. Since the position comparator remains 
enabled the increase in the the position error is due solely to direct interference of the “protons” with the 
star, since the “proton” signals merge with the star signal to form an erroneous star position centroid. 
Despite the relatively large position errors and the numerous dropouts of “star’ data, the tracker never lost 
the “star” to the extent that it re-entered the acquisition mode, see View (b) of Figure 9. 

Figure 10 presents the results from the same SS test scene but with both the intensity compare and position 
compare discrimination of the tracker disabled during the track mode. The standard deviation of the reported 
“star” position increased to 145.3 and 173.7 arc-seconds in row and column positions, respectively. These 
large errors are due not only to interference with the “star” by the “protons” but, in addition, “proton” gen- 
erated centroid positions are mistakenly reported as the “star” position. This case of mistaken identity 
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Table 1. Average Acquisition Times in “Proton” Rich Environment 
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Figure 8. Tracker Performance With Proton Flux Design Features, (a) “Star” Position Variations During 
10 Minute Test; and (b) Row/Col/Intensity Values Plotted Over the Test Time. 
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Figure 10. 


Effect of “Proton” Events on Tracker Performance Without Proton-Flux Design Features, 
(a) “Star” Position Variation Over 30 Minutes Caused By “Proton” Events; and (b) Part of 30 
Minute Test. 

occurs if the proton signal occurs within the Track Window and precedes the star position during the raster 
scan of the window. The dissymmetry of the “star” position variations shown in view (a) of Figure 10 is 
explained by this phenomena. View (b) shows a case in which the errors have become so large that the 
“star” is lost and the tracker must re-enter the acquisition mode to re-acquire the star . 

Figure 11 evaluates the same data from an acquisition perspective. The software of the tracker with t the po_ 
sitfon and magnitude comparators disabled was also modified so that during acquismon the first signal with 
the characterises of a valid star (i.e., proper size and magnitude) would cause a transition to track^ (The 
multiple frames and tests implemented in the design for use in high-proton flux environments were disabled.) 
In the proton-rich environment, as simulated here, a large number of false star acquisitions were cause y 
proton signals. The actual “star” was acquired only nine times in 30 minutes of tracker operation, and was 

tracker for only six percent of the time. 
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Figure 11. Without the Proton Flux Design 
Features, the Tracker Will Usually 
Acquire and Attempt to Track 
Proton Events. 


SUMMARY 


Radiation-induced proton events can result in anomalous operation of solid-state star trackers, specifically: 

• Erroneous acquisition of proton events and/or failure to acquire valid stars 

• Loss of track of valid stars 

• Incorrect position and intensity data. 


Through incorporation of hardware-efficient processing algorithms, HDOS has completed delivery of two flight 
trackers for the TOPEX/Poseidon mission which can operate effectively in a proton-rich environment. Using a 
scene simulator to produce effects similar to those caused by protons, tests validating the performance gains 
achieved have been completed on both units. For an environment that produces 150 false multi-pixel events 
at the detector, the following results were obtained: 

• No acquisition of false stars (proton events) 

• Reliable acquisition of valid stars 

• No loss of tracking a valid star 

Identification of corrupted data for the host, caused by proton impact upon valid star pixel groups. 


The algorithms incorporated into the tracker firmware can be tailored to unique user mission requirements 
The scene simulation techniques developed provide a powerful tool for validating performance for rather 
unique and complex test conditions. 
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